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Abstract 

We present a novel approach to Gaussian Berezin correlation func- 
tions. A formula well known in the literature expresses these quantities 
in terms of submatrices of the inverse matrix appearing in the Gaussian 
action. By using a recently proposed method to calculate Berezin inte- 
grals as an expectation of suitable functionals of Poisson processes, we 
obtain an alternative formula which allows one to skip the calculation 
of the inverse of the matrix. This formula, previously derived using dif- 
ferent approaches (in particular by means of the Jacobi identity for the 
compound matrices), has computational advantages which grow rapidly 
with the dimension of the Grassmann algebra and the order of correlation. 
By using this alternative formula, we establish a mapping between two 
fermionic systems, not necessarily Gaussian, with short and long range 
interaction, respectively. 

1 Introduction 

The concept and the use of anticommuting variables originates both in the con- 
test of the functional integral approach to the quantization of fermionic systems 
[1, 2] and in several combinatorial problems (e.g. to represent the partition 
function of the planar Ising model) [3, 4]. The anticommuting character implies 
that these variables belong to a Grassmann algebra. 

A Grassmann algebra Q g of g generators {£,i,£,2, ■ ■ ■ ,£, g }, is defined by the 
identity l 9 and requiring the generators obey the following anticommutation 
relations: 

= 6& + &6 = ° '•./ • {i- •••//!• (i) 

Elements of Q g of the form ChS,i 2 ■ ■ ■ £ik with 1 < ii < *2 < ■ • • < *fc < <? and 
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k < g, are called monomials. Each elements of Q g , F(£), can be written as a 

unique complex linear combination of the 2 9 (independent) monomials: 

g 

F«) = f i g E 

fe=l l<ii<i2<---<ifc<g 

where fi lt i 2 ...,i k G C. In this algebra, derivation and integration (the Berezin 
integral) are properly defined. The reader is referred to [1, 5, 6] for details. Let 
us introduce the symbols d£i, . . . , d^ g satisfying the following formal relations 

{d£M = {dZ i ,S j } = 0. (2) 
Let us define the elementary Berezin integrals 



/ 



B pB 

d£i = 0, and for j =/= i / </i, = 



iterating this relations we arrive to 

so, from linearity we have the following general definition of Berezin integral 



(3) 



The derivation is defined as 

+(-l) n - 1 WU^ 2 "-^»-i- ( 4 ) 

Let now C7 2 g be a Grassmann algebra over C, with 2g generators, which we 
label as g barred and g unbarred: £1,62, • • • , £ 9 ; £1, £2, • • • , £»• Let Z 2 = {0, 1} 
and (Z 2 ) X9 = Z 2 ® Z 2 ® ••■ ® Z 2 ( tnc direct product of Z^ 5-times). Let us 
define the normalized Gaussian Berezin integrals: 



/ 



where A is any invertible g x g matrix, ^ = (£ 1; £ 1; £ fl , £ fl ), <0 S = £™ 1 ^" 1 . . . £g 9 ^ ! 
with 

s = (ni,ni, ...,n g ,n g ) G (Z 2 *) 2s (6) 



and the definition of the Gaussian measure dfi is implicit in (5). Let k be the 
number of labels i such that fn = 1 and k the number of labels j such that 
rij = 1 . Due to the parity of the Gaussian function exp(— f, A£), it is easy to 
see that the integral (5) vanishes if k ^ k. 

When k — k, the integral (5) can be evaluated by a well known formula in 
terms of the determinant of a submatrix of the inverse matrix C = A -1 (see 



where C s is the matrix obtained from C erasing the rows with label j if nj = 
and the columns with label i if fii = 0. If s = (0, . . . , 0), we define C s — 1. 
Equation (7) is then expressed in terms of a k x k submatrix of the inverse of 



Starting from a recently developed approach to calculate Berezin integrals 
[7], we will derive (see Eq. (12)) a formula for (5) in terms of the determinant of 
a (g — k) x (g — k) sub matrix of A. When compared with the (7), this formula 
gives rise to an algebraic identity, which is nothing else than Jacobi's theorem 
relating the minors of a matrix to the minors of the inverse of the same matrix, 
see e. g. [8, 9], see also [10]. Although this identity and the Gaussian Berezin 
correlation functions can be derived algebraically in a simpler manner, even not 
using the Jacobi's theorem [11, 12, 13], a probabilistic derivation provides an 
important example of how the stochastic approach for fermions models can be 
applied to obtain exact results. 

It is worth considering that even if this identity is well known among math- 
ematicians, it seems that it has been rarely explored in Quantum Field Theory 
[11, 12, 10]. In section 5 we will explain some interesting aspects and physi- 
cal applications of this identity. In particular we will find a mapping between 
two fermionic systems. According to this mapping the Grassmannian partition 
function of a short range system having A as matrix for the free (Gaussian) 
part of the action, equals the Grassmannian partition function of a long range 
system having A" 1 as matrix for the free part. 



e.g. [2]): 




(7) 



A. 
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In the remainder we will give the proof providing in full detail all the calcu- 
lations as they may look unfamiliar to many readers. 

2 Gaussian Berezin integrals 

Given s = (m, h\, n g ,n g ) G (Z 2 *) 29 , let 

s* = (1 - ni, 1 - ni, 1 - fig, 1 - rig). (8) 

Probabilistically, we will show that besides 
/• B - 

/ exp(-£, A€)il> s d£id£ 1 ...d£ g dt g = detC s det A (9) 

holds also 

y exp(-€,^)^ s deidei»-de fl de fl = det^.-(-l) W(S) - 1 ( 10 ) 
The factor(— l) w (s) determines a global sign with 

W{s) = (m + ni) + 2(n 2 + n 2 ) + • • • + g(n g + n g ) (11) 

From the above equations we get the following expression for the normalized 
Gaussian Berezin integrals 

and by using Eq. (7) the Jacobi's identity follows 2 

3 Berezin integrals as averages over Poisson pro- 
cesses 

Let Z 2 = { — 1,1} and Z% 9 = Z 2 ® Z 2 ® . . . <g> Z 2 (the direct product of Z 2 g-times) . 
Let e and er G 7L^ 9 and e ff = (1, . . . , 1) G 9 . Let {-/Vf} e ^e s be a family of 

According to the definition given for C s in introduction, A s * is the (g — k) X (g — k) 
submatrix of A obtained from A erasing the rows with label i if fii = 1 and the columns with 
label j if nj = 1. If s = (1, . . . , 1) we define A s * = 1. 

2 It is simple to observe that (ni + ni) + 2(n2 + «2) + ■ ■ • + g{n g + n 3 ) = ii + ji + «2 + 
J2 + • • ■ ifc + ifc , where ii . . . jj, and ji . . . j'fc are the the sets of indices such that = 1 and 
rij„ = 1 respectively, for v = 1 . . . k. In this form the sign in the Jacobi's relation is normally 
presented in literature. 
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2 9 — 1 left-continuous independent Poisson processes of pure growth with unit 
parameter. These Processes are then characterized by the probabilities 

P(Nf +At -Nf = k) = ^f- exp (-At) (14) 

Following [7] 3 , given any "action" S(£) E Q g , we can calculate the Berezin 
integrals by an expectation of a suitable functional of Poisson processes using 
the following formula: 

/V •■• ^»exp(-5(€))dei...de fl = A fl (-(-l)")exp(|r|- S (e fl )) 
yE ( Jll^-(-l) Nl i +n ' 

x exp j]T ^ log(-*(e)C fl (e, -(-1)^ + n ))dNf^ ,(15) 



where: 



i=i 



A g (er) = n( V^ + ^-i-.-^-i), (16) 



s 



^ *) = n + • • • • • • ^-i j . 

(-i) n = ((-ir,...,(-i)"')ez 2 x9 , (18) 



s(e) = A(e) / + €) 5(€)dei •• where 

1/1 J,-, 1 1 (l+«g) 

^(1 + 6) = ^ _£ g2 ); (ig) 



r = {eeZ 2 X5 :e / e fl , s(e)C g (e, <x) ^ 0}; |r| = the cardinality of T, (20) 

N t = (Ni,...,N?); Nl = Y,\^-ei)Nf. (21) 

eer 

Finally the stochastic integrals that appear in (15) arc ordinary Stieltjes inte- 
grals: 

/ f(8,N a )dN a =J2f(8h,N. h ), (22) 

J\o,t) t~i t 



Sh<t 



3 In [14] it is shown in a simple manner, how to derive the Poissonian probabilistic repre- 
sentation of the time evolution operator for any operator described by a finite matrix. This 
derivation can be used instead of that appearing in the fifth chapter of [7]. 
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Sh being the jump times of the process. We refer the reader to [7, 14] for proofs 
and details. Here we want to emphasize that: only the trajectories such that 
Nf £ {0, l},Ve G T can contribute to the expectation in (15). We will indicate 
such a set as the set "zero — one" . This happens due to the structure of 
the coefficients C g present in the argument of the logarithm which occurs in Eq. 
(15). In fact, let us consider the expectation of the exponential occuring in (15), 
if we look at the definitions of the coefficients C g , neglecting their moduli, we can 
write this expectation as a sum of terms corresponding to the following families 
of events with respective contributions and probabilities (the latter are obtained 
by integrating over the jumps times the infinitesimal probabilities (14)): 

• event in which no process jumps. This event gives 1 with probability 
exp(-|r|); 

• events in which only one process, say with a generic label e 1 , makes a 
jump. Each of these events will give the value — s(e 1 ) J[ l . e i = _ 1 ' 
with probability f Q ds\ exp(— |r|) = exp(— |r|); 

• events in which two and only two processes with different labels, say e 1 
and e 2 , jump at times s\ < s 2 respectively. Each of these events will give 

the value (-«(^))n i!ef= _ 1 * (--(^IL^-i ^T^^ 

with probability s ^ <s ^dsids2 exp(— |r|) 

• and so on . . . 

• events in which one single process, say with label e 1 , jumps two times. 
Each of these events gives: -sfe 1 ) Ui-^=-i 1+{ ~ ir ' x (-^(e 1 )) Ui-.e^-i 
with probability J Q s ^ <s ^dsids2 exp(— |T|) 

• and similarly for processes with more than two jumps. 

It is then clear that the last type of contributions gives always zero. Furthermore 
note that events in which two or more processes jump at the same time, have 
probability zero, even if C g ^ 0. 
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So, given S(£) and then s(e), one will have to determine what is the set of the 
effective trajectories, in general a smaller subset of the above set of trajectories 
zero — one. In the next section we will investigate the Gaussian case (S bilinear 
in the £) and we will see that such a set is "very small". The search for the 
effective trajectories and their phases will allow us to automatically construct 
the matrix A s * and, in the final expression, to recover its determinant up to a 
global sign. 

4 Proof of equation (12) 

The application of Eq. (15) will give us the Gaussian integral not normalized, 
let be I (the numerator of (5)). Let us consider (15) with g barred variables and 
g unbarred variables. We will have to calculate various parts. In section 4.1 we 
will calculate s(e) and then we will have T for our bilinear case, next we will 
calculate partially Ci g - This result will be used in section 4.2 for calculating the 
effective trajectories and in section 4.3 for evaluating the phase of (15) which 
will provide us the final formula. We will consider the integral with respect to 
d£i...d£ g d(;i...d£g (instead of d£id£i...d£gd£g). Therefore, we will evaluate the 
following integral 

1 = j B .,«:...</^r • • ■ ■ ■ ■ Q 9 exp(-|, M), (23) 

which is related to I by a simple correcting sign, / = I(— 1) Q which we will 
calculate in section 4.3. We will see that Q depends on g and s. With the 
chosen ordering of the generators of the Grassmann algebra, e and s should be 
changed into 

e= (ei,...e g ;ei,...,e fl ), (24) 

8 = (ni,...,n g ;ni,...,n g ). (25) 

However, for simplicity, we will continue to use, if not ambiguous, the initial 
notations e and s also for the new arrangement of the Grassmann generators. 
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At the end we will return to the true integral corresponding to the original 
arrangement of the Grassmann generators by using I = I(— 

4.1 Calculation of the factors: s(e), A 29 (e) and C 2g 

Let us start by calculating s(e) defined by Eq. (19): 

s(e) = A 2g (e) ' . . . ti^k^ ■ • • ZZ^fc M) 

= A 2g (e) f ( /:,..,/C ; ^ l ...< ; ^ ' . . .# (1+f9) ^ (1+ei) • ..a (, " ! ' ; ^>^ 



= A 2g ( e ) 5; / B dei...de fl dei...de fl ^ (1+El) . . . . . . 

ij J 

xrf (1+ei) ...e/ (1+£j)+1 ...a (1+£9) , (26) 

where (— l) Pij accounts for the exchanges needed to bring & in front of £?^ 1+e ^ 
and £j in front of £?( 1+ej ). We see, that integral (26) vanishes unless e has the 
form 

eij = (ei, ...E g ;ei...,e g ), with E l = ej = -1 and ej = ej = 1 for Z ^ (27) 

Furthermore it follows from this that 

s(e 2s ) = 0, where e 2ff = (1, 1; 1...., 1). 

The set T is then characterized. Its elements are in the set of the 2g — pie €ij 
such that Aij ^ 0. In this way we can arrange the set of Poisson processes in a 
g x g matrix: 

Af = ((Nfj), where N? = . (28) 

Even if not explicitly indicated, all the sums and products will have to be 
considered in the range such that ^ 0. 
From Eq. (26) we have 

s(e ij ) = A 2g (e ij )(-l)^A ij . (29) 

In appendix B it is easily proved that one has 

s(e ij ) = A ij . (30) 



As regards the calculation of Ci g , according to our choice of the ordering of the 
Grassmann generators, we must consider (see Eq. (21)) N s = (N S ;N S ) and 
n = (n; n), where 



eer 



N l s = ^hl-e l )Nf. 



(31) 



(32) 



eer 



So to calculate Ci g one has to evaluate 

C 2fl (e«,-(-l)^-+") = C 2 ^ 

In appendix C it is easily proved that one has 

c ,( £ „,-(-D"--) - 1+( -f*"' 1+( -f"V ir'+'+- 

x exp{i7r[AT s l + ... + N 9 a + N$ + ... + N*' 1 + n t + ... + n g + n x + ... + rij-i]}, 



where 



1=1 



(33) 



i=i 



Therefore, as far as the exponential factor in (15) is concerned, a given trajectory 
will give the following contribution 

exp / \og(-A l3 )dN^ exp / log[- 1 ' 
i:j U[o,i) J U[o,i) 



x exp 



2 



/ »7r[ ff + i + j + 1 + (NT + rim) + £ (N? + n m )]dN? 

-"'I . 1 ) m=i m=l 

where the range of the product is over the couples ij with ^ 0. 
4.2 Effective trajectories 

In section 3, we have seen that the set of the effective trajectories is a subset 
of set zero — one (i.e. the set in which one single process can make at most 
one jump during the time interval [0, 1)) and we have anticipated that in the 



]dN? 
(34) 



Gaussian case such a set is small. Let us concentrate on the second factor of 
the product (34). Let us define: 

G 2 -G, = — ^ ±L_J (35) 

and remember that the processes are left continuous 4 . Let us note from (33) 
that iVj and N£ are the sums of the elements of the ith row and the j'th column 
of the matrix A/" respectively. 

Let us now consider the trajectories in which in every row and in every 
column of Af there is at most just one process which jumps (clearly, one process 
which jumps in some row, jumps in some corresponding column too). With such 
events, if for example N l J is a process which jumps at time s*- 7 , then according 
to Eq. (35) we have: 

• Lrj|s=s i J — ^ 2 ' 

Since we consider the factor corresponding to the link (ij) and carry out the 
stochastic integral in dN l J , we see that in order it does not give explog(O), it 
must be fn = rij = 0. 

Let us now consider one event in which there is at least one row in which at 
least two processes, N l J and N lk , jump, at times s y and s lk respectively. Let 
us consider the (ij) factor in Eq. (34). We have to integrate in dN l J , i.e. we 
have to evaluate the integrand at time s u . Let fii = 1. We must consider two 
possibilities: 

1. if s y < s lk we have: 

N l s ij = from which follows Gi = 

2. if s y > s lk we have: 

N l gij = 1 from which follows Gj = 1. 

Therefore the stochastic integral can be finite in the the case 2, but as we will 
integrate in dN lk , we will have to calculate the integrand at time s lk then ob- 
taining Gi — 0. Finally Let us note that the event s 4 - 7 = s lk has zero probability. 

4 This implies that if s = s <J is the jump time of the first jump for the process Nl J , then 
Ni J = for s < s ij and N' s j = 1 for s > s ij . 



10 



Let fii = 0. Repeating a similar argument we will arrive at the same conclu- 
sions. Of course, the same discussion can be done for the columns. 

We have proved that the effective trajectories belong to the set in which 
in any row and in any column of the matrix Af there is at most one jumping 
process. Let us call such a set: 

Ci = {set of trajectories in which for any row and for any column 

of the matrix Af there is at most one jumping process}. (37) 

Furthermore we have proved that as fn = 1 and/or rij = 1, then it must be 
N\ 3 = 0. Let us call such set: 

C-2 = {set of the trajectories in which if fii and/or rij = 1 

then N{° = 0}. (38) 

If now we consider (34), (regardless of the phase), imposing these restrictions by 
using the product of two characteristic functions, we can eliminate the stochastic 
integral writing (34) as: 

n(-^)^-x(Ci)x(C 2 )<% (39) 

where with .S^ we have indicated the phase (i. e. the last of the three factors 
appearing in (34)). Observing Eq. (15), we see that to complete the search of 
the effective trajectories, we must multiply the above product by the following 
factor 

A i - f-i)^+™' i - (-i) N l+ n i 

F - F = Ii L ^ —T • ( 4 °) 

This factor takes the values or 1 and its presence, using the fact that in Eq. 
(39) the function x{C\) occurs, involves a third characteristic function x{Cs) 
which accounts for the fact that if for some to, n m = 0, then there must be 
one and only one j such that N™ 3 = 1. Analogously, if for some m, n m = 0, 
then there must be one and only one k such that N± m = 1 . If now we use also 
the function x(^2), then for any given s (6) we find that the set of the effective 
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trajectories C is given by: 

C = {set of trajectories such that: (41) 
if fii = 1, any process on the row i of J\f does not jump; 
if rij — 1, any process on the column j of does not jump; 
if fii — and rij — 0, there is one and only one process jumping to 1 
on the ith row and on the jth column of A/"}. 

Clearly if k ^ k then C reduces to the null set and the integral (23) gives 0. 
From now on, we will always suppose the k = k case. 

From its definition we see that C has (g— k)\ elements so that we can represent 
it in terms of suitable effective matrices Af in the following way. First we must 
erase all the rows with label i if n» = 1 and the columns with label j if nj = 1. 
Then we must consider all the ways to fill the remaining (g — k) x (g — k) sub- 
matrix with (g — k) times "1" in an array such that two of them are never in 
the same row or in the same column, like for a determinant. 

Let us consider as an example the case g — 5 and k = 2 with n\ = n 2 = 
ns = fii — fii = fiz = and n 4 = n 5 = n 4 = fi 5 = 1. In this case C has 6 
elements represented by 6 matrices A/", like e. g. (the symbol "x" means that 
the corresponding process never jump) 



( 


1 





X 


X 


\ 


1 








X 


X 










1 


X 


X 




X 


X 


X 


X 


X 




V >< 


X 


X 


X 


X 


/ 



Up to now we have not spoken about the event in which no process jumps 
(with probability e~' r '). It is immediate to observe that such an event is effective 
only if rii = hi = 1 for any i from 1 to g and that the integral (23) gives simply 
1. So, from now on, we will ignore this trivial case. 
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4.3 Phases 



Let us face now the problem of finding a closed expression for the Sij (see Eq. 
(39)): 



S tj = exp 



/ in[g + i + j + l + J2(N? + n m ) + £ (N, 

■ J M m =r m=l 



(42) 



Then we will need to consider the product on the couples (ij) for every (ij) 
such that N\ 3 — 1. Let us factor the rhs of the above equation in the following 
form: 



x exp 



m[g + i + j + l]Nl 3 

3-1 



Sij = exp 

/ (NT + n m ) + f^ (NT + 4)]^ 

^[o,i) m= i 



(43) 



Let s lJ , s lm , s ml be generic jump times with their corresponding labels, then 
from (33) we obtain: 



= exp 



in[g + i + j + l]Ni j 



x (-1) F ^ 



(44) 



where 



g a 



Fij = [J2J2 e ^ 13 - sml ) N r l + 2 e ^ l ° - slm ) N i m \ N i + 

m—i l — l m—1 1=1 

g j — 1 



m— 1 



which, taking account the definition of C\, becomes: 



m=i+l /— 1 m— 1 (=1 

g j — 1 



Now in order to calculate the expectation, we must integrate over the jump 
times and then to sum over all the trajectories the following product: 



n(-^)x(C)exp 



in[g + i + j + l]Ni j 



x (-1) F «. 



(47) 
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The integration must be done using the probability of the Poisson processes and 
it is easy to see that 

dP = e-^ Y[ dsiJ - ( 48 ) 

y:Af =1 

We should integrate (47) with the measure dP and then sum over all the possible 
events belonging to the set C; we should then calculate the following integrals: 



L 



Vl r l 

o 



j] {-l) Hi 'ds ij , (49) 

i :Nl 3 =l 

where 

H i:j = [ E E °( sij - s " d ) N r l + E E 9 ( sij - slm ) N i m }- 

m=i+l 1=1 m—1 l—l 

In appendix D it is proved that 

= [ E E c^r' + E E « m ] . (50) 

m=i+l 1=1 m=l 1=1 

where the variables with four labels 9 l ^ k satisfy the following properties: 

(51) 

Now we have the following formula for the integral I (the numerator of (5)) 

/ = (-i)«A 2fl (-(-i 5f )) yi •• e ^ n (-^)(-i) 9+i 

x (_i)EL i+ i * m +Ei=, «m , (52) 

where s = (n, n) and (— 1) Q is the sign necessary to have I from 7. 

Now we will face the problem to obtain a closed expression for each of the 
factors appearing in the above equation for I. 

4.3.1 Correcting sign (1)^ 

To calculate (— 1)^ we must consider the exchanges needed to restore the mono- 
mial . . . fg 9 ^™ 1 . . . £,g B and the ordering of integrations in the original order. 
The first produces a sign equals to: 

(— l)™i(™i+™2 + ---+™s)+ n 2(™2 + ---+ng) + ---+n g n s 
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the second: 

(_l)fl(fl+l)/2_ 

Therefore 

(_1)<2 = (_ J^ni(ni+n2+---+n 9 )+n2(n2+---+S<,) + ---+n a S 9 +9(s+l)/2 ^gg-j 

4.3.2 A 2s (-(-l 3 )) 

For A 2g (-(-l 5 )) we have: 



-(-if 1 ...(-l)™'- 1 



1=1 

X Y[( V, + TT^(- 1 )™ 1 • • • (-l)" 9 (-l) ni • • • 

1=1 ^ 



2 2 
from which one easily has: 



A 2s (-(-i ? ))=n ex p 



ni H h n;_i + / - 1 



PJexp^(i-7r)nZ 



i=i 



ni H h n g + m H h n;_i + <? + / — ! 



So we have: 



A 2ff (-(-l 5 )) = 



(54) 



where 



i? = n 2 {ni + 1) + n 3 (ni + n 2 + 2) H hn g (n 1 H h n s _i + # - 1) 

+ni (fc + .9) + n 2 (fc + .g + ni + 1) + . . . 
+n g (k + .9 + ni H h + .9 - 1), 



i.e.: 



R = n 2 {n 1 + 1) + n 3 (hi + n 2 + 2) H hn g (niH hn r i+g-l) 

+n 2 (ni + 1) H h n 9 (m H h + .9 - 1) 

+(fc + s)(ni + ---+n g ). 
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4.3.3 Constant factors 

Due to the structure of the set C, it is easy to observe that the following relations 
hold: 

x (c) 11 (-i)SU.*»+Ei,=i% = (_i)S 

ij: JV«=1 

with 

S = (l-ni)(n 2 H h n s ) + (1 - n 2 )(n 3 H h ra s ) + . . . (1 - ra s _i)(n 9 ) + 

(1 - n 2 )(n 1 ) + (1 - n 3 )(m +n 2 ) . . . (1 - n g )(m + . . .n g _i) (55) 

and 

X (c) n (-ir +j = (-if, 

ij: N^=l 

with 

T = ^T(l-n^ + ]T(l-n,)j. (56) 
Bringing out of the sum the constant factor (— 1) 9+1 (— 1), we arrive at 

I = (-l)Q^ + T +g(9 - k) £ £ x(c) -q (A^) (-1)^.(57) 

Now we must calculate the parity of the sum Q + R + S + T + g(g — k) . We will 
indicate the parity of an integer n with: 

V[n] = Parity [n]. 

Using the fact that J2i ™i = J2i n * = it is not difficult to get the following 
relations: 

V [R + S + T] = 7>[(m + nt) + (n 2 + n 2 )3 + • • • + (n g + n g )(2g - 1) 

+n 1 (n 2 H hS 9 )+ n 2 (ni + n 3 H hn g ) + ... n s (rii H h "g-i) 

+(n 2 H hn s ) + (n 3 H 1- n g ) H h (n s ) 

+(ni) + (ni + n 2 ) H h {n\ H h n s _i) + k(g - k)], 
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from which one has: 



V[R + S + T] = V[{n 1 +n 1 ) + {n 2 + n 2 )2 + --- + {n g +n g )g 



+k(k 



l) + k(g-k) + k(g-l)], 



so, using Eq. (53) we obtain: 



V[Q + R + S + T + g(g 



k)} = 



P[ni(ni + n 2 + ... + n g ) + n 2 (n 2 + ... + n g ) + ... + n g n g 




In the next subsection we will extract another constant factor coming from the 
phases (—l) Hij contained in the product (57). 

4.3.4 Permutations - completion of the proof 

Let us now investigate the full meaning of the H^. Let be h = g — k. Let X 
and J be the sets (of cardinality h) of the labels i and j satisfying fii — and 
rij = respectively. We will call such sets "active" . First, let us suppose that 
X = J = {1, . . . h}. In terms of the matrix A, we are looking at the submatrix 
of A obtained considering the first h rows and the first h columns of A. 

In section 4.2 we have seen that each event of the set C is a proper realization 
of the matrix Af. Choosing the rows as reference, we can obtain the events of C 
considering for each of the first h row indices i, a permutation 7r = (m, . . . , -kh) 
such that N™ i = 1. In other words, for any row i £ X, ni tells us which is the 
corresponding (unique) active column. With these definitions we can carry out 
the summation over the events by summing over the permutations: 



where II is the set of permutations of h elements. Let us define also for each of 
the first h column indices j, the corresponding permutation u> = (wi, . . . ,Uh) 
such that N^ 3J =1. If tt and lo are referred to the same event, it is easy to see 
that one has u) Vj = j for any j £ J. 



E *(c)/(-) = E /(•). 



.,JVf=0,l n II 
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Given an effective trajectory, if ir is the corresponding permutation, one has 

h h TTi — 1 

Y Hi i = Y [ Y 6 ™m + Y 6 ™i™\- 
ij:N^ — l i=l m=i+l m=l 

Let us now expand the r.h.s. of the above expression in the following way: 

[elii + + ■■■ + C] + [Ci +«& + ■■■ + Ct.-J 

+ + • • • + CJ + [Ci + Ca + ' ' ' + C 2 2 - 2 )(--2)] 
+ 

+ 

+ 

Let us consider the first row of the above expression. If 7Ti > 7r 2 , then in the first 
row, besides 6> 2 ^, 0^J^ n2 occurs too and w V2 = 2. On the other hand if m < 7r 2 
(the equal sign is not possible since 7r is a permutation) the former pairing for 
9^1 is not possible in the same row, therefore the only possible pairing is with 
a different row where the upper labels for the 9's are all different from lwi. So 
in this case the pairing will be with 6\l\ . Then, using the relations (51) one has 
the following dichotomy: 7Ti > 7r 2 — > 1 and 7Ti < 7r 2 — ► 0. Repeating the same 
argument for any row of the above expression one arrives at: 

J2 H i:j = 9(n 2 - tti) + 9(n 3 - tti) + • • • + 0(ir h - m) 

+0(n 3 - tt 2 ) + • • • + 6(n h - tt 2 ) . . . 

+ 

+ 

+.... 

+9{n h - 7T h _i). (60) 

The r.h.s. of (60) is just the number of exchanges needed to bring the per- 
mutation 7r to the inverted fundamental permutation i.e. to the permutation 
(M- 1, •••,!)■ 
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If now we want to express the sign of the permutation in more conventional 
terms, i.e. in terms of the number of exchanges needed to bring the permutation 
to the fundamental one (1, 2, . . . , h), it is enough to consider the relation between 
such number of exchanges and the former. If we define Hij in the following way: 

J2 = [i_e( 7r2 _ 7ri )] + [i_e( 7r3 _ 7ri )]+...+ [i_e( 7rfc _ 7ri )] 

ij:Ni j =l 

+ [1 - 6(lT 3 - 7T 2 )] + • • • + [1 - 0(TT h - 7T 2 )] . . . 

+ 

+ 

+ .... 

+ [l-e(TT h -W h -l)], 

we see that the above r.h.s. now counts the number of exchanges needed to 
bring the permutation to the fundamental one; on the other hand we have: 

£ H ij= £ Hij - h ^A. (61 ) 

ij:Nl j =l ij:N{ j =\ 

The term MftiziJ j s the last of the constant factors just mentioned before. It 
must be summed to all the others constant factors, i.e. to the r.h.s. of Eq. (58). 
So under our hypothesis I = J = {1, . . . h}, we have the following final global 
phase: 

(_ l ^[Q+B.+ S+T+g(g-k)+!i<± T ll} = ^_^W (g 2 ) 

where Wo is given from 

Wo = ni(ni + n 2 + ... + n g ) + n 2 {n 2 + ■■■ + n g ) + ... + n g n g 

k(k + 1) 

+(ni + hx) + (n 2 + n 2 )2 + • • • + (n g + n g )g + — - 

On the other hand if 1 = J = {1, . . . h} it is easy to see that (~l) w ° = 1. 
Now we can return to the general case in which the active rows and columns 
are no more restricted to be the first h. For this aim it is enough to consider 
the numerator of the (5) and to make a suitable permutation which brings 

19 



d^idii-.d^gdig in d£ Pl d£ qi ...d£ Pg dt; qg with n pi = n qi = . . . n Ph = n qh = 0. So we 
have a again the case I = J = {1, . . . h} up to a sign related to this permutation 
which is (-l)E^=iP-+^ = (_l)ES=i^+^ = (-1)^ (see the foonote n.2). 
Then for the Gaussian, not normalized, Berezin integral / we have found: 

i = (-i) w £(-i) ff n^ 

where now 7r = (ir\, . . . ,-Kh) with ij e J for any iel and (— l) 71 " is the sign 
of the permutation 7r, i. e. the sign determined from the parity of the number 
of exchanges needed to bring 7r to . . . ,jh) with j\ < ■ ■ ■ < jh- So up to the 
global phase (—l) w , I is just the determinant of the (g — k)x(g-k) submatrix 
defined in section 2. Then we have found: 

7-dctA s .(-l) w ( s ) 

5 Algebraic, geometrical and physical aspects 

Equations (12) and (13) provide interesting interpretations and applications in 
several contests. 

5.1 Algebraic 

If we choose s = (1, 1, . . . , 1, 1), Eq. (13) reduces to the trivial identity 

detC=^, (63) 

which tells us, as it is well known, to calculate the determinant of the inverse 
of A, it is not necessary to evaluate the inverse of A. Eq. (13) generalizes this 
property to any submatrix of C, C s , since to obtain detC s we need just the 
determinants of A and A s ., with A s . submatrix of A. However it is important 
to note that in general C s (— l) w is not the inverse of A(A S * where with 
B s we indicate the g x g extension of JB S , i.e. the matrix representing the 
operator which is the direct sum of the operator in k dimensions, represented 
by the matrix £? s , plus the identity operator in g — k dimensions. 
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Actually we can easily prove algebraically Eq. (13) in cases in which A a * is a 
submatrix of A with contiguous labels, in the following way. Let us suppose A s * 
is the submatrix obtained erasing the last h rows and h columns of A, where 
h = g — k. Let us decompose the matrix A in four blocks: An; A 12 ; A 21 and 
A22, with dimensions kxk; kxh; hxk and hxh respectively. Let us make the 
same for the matrix C. In particular we have A 22 =A a * and C\\~C a . With 
such a decomposition it is easy to multiply any two matrices by considering 
the blocks like elements of matrices of dimensions 2x2. Now let us define the 
matrix B=A S * + D, where D has £>22=^4i2 and the others three blocks zero. 
Hence one has det £?=det A s *=det A s * . On the other hand for E = CB one 
has £ ; n=Cn; E 12 =0; E 21 =C 21 and E 22 =l. So, since det E = detC s , by 
using 

dctA s . , 
lMA- = dct ^ 

we get the identity (13) (it is simple to see that in this case W is even). 

In the general case we will have to face the complexity of the indices and 
the proof of Eq. (13) can be obtained as a corollary of the Binct - Cauchy's 
theorem on the compound matrices [8, 9]. 

5.2 Geometrical 

Up to a sign, the determinant of a g x g matrix represents the volume of a prism 
in a g dimensional Euclidean space. The volume is that spanned by g vectors 
whose g components are the g rows of the matrix. From Eq. (13) we have 

' detA ' = ^— (64) 
|detA a .| |detC s | v ; 

Now, the l.h.s. of this equation is the ratio of two volumes: one is that of a 
given g dimensional prism, while the other is that of a g — k dimensional prism 
obtained by spanning the space using the g — k vectors whose components are 
given by the g — k rows of A s * , i.e. it is a g — k dimensional section of the first 
prism. So this ratio is nothing else than the relative "height" with respect to 
this section (here with "height" we mean a generalized height as in general it 
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will be itself a lower dimensional volume). Hence Eq. (64) tell us that | det C s \ 
is the inverse of this height. 




H = 1/| det(C a )| 



S= |dct(A a *)| 



5.3 Physical 

5.3.1 Gaussian averages 

Equations (7), (12) and (13) show that the problem of evaluating Gaussian 
Berezin integrals can be addressed as follows: 

• if A^ 1 is known, then we will use Eq. (7), 

• if A is known, then we will use Eq. (12), 

• if both A and A^ 1 are known, we will use Eq. (7) if k < g/2 and Eq. 
(12) if k > g/2. 

The last point is obvious. Since the calculation of a determinant of a d x d 
matrix leads to a sum of d\ elements with alternate signs, the advantage in 
the choice of one or another representation may be factorial with respect to 
the dimensions as k is enough different from g/2. The others two points, first 
remarked in [11], are also obvious. In fact let us suppose we are, for example, in 
the second case and we are interested to the not normalized Gaussian integrals 
(that is, we don't need to evaluate det A for both the two representations). In 
order to calculate C s from A we must calculate k 2 minors of A, each of them in 
general involving a determinant of a (g — 1) x (g — 1) submatrix of A. Therefore 
we should to consider k 2 [(g — 1)1] terms in order to calculate detC s , while to 
calculate det A s * we need to consider (g — k)\ terms. So, in general, for any 
k > 1, knowing A, the use of Eq. (12) provides the following gain (calculated 
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as ratio between the number of terms present in the two representations) 

GAIN = k 2 (g - l)(g - 2) . . . (g - k + 1), (65) 

while for k = 1 the GAIN is 1. 

For normalized Gaussian integrals, actually the gain can be less as the calcu- 
lation of det A could implicitly provides some of the above k 2 minors. However 
in general the gain will be of the same order. 

5.3.2 Mapping between two fermionic systems 

Besides the above practical advantage, Eq. (13) provides an interesting theo- 
retical insight. 

In terms of Gaussian Berezin averages Eq. (13) may be read as: 

(_l)W(S) 

< r > A =< v > A - 4cb— (66) 

< ' > A means a Gaussian average being A the matrix of the Gaussian action. 
Let be given a fermionic system with the following action 

= \(l M) (€,€), (67) 

where V(£, £) is an arbitrary element of Q 2g which we can expand in monomials 
as 

SG(Z 2 *) 2s 

with Vs € C. The system is then characterized by the following Grassmannian 
partition function 

r-B 



Z = J &q>[-S(t,t)]dZ 1 d$ 1 ...dZ g d€ g . 



(69) 



Let s* = (mi,mi;...;m 3 ,m s ). Let f be a linear operator on Q 2g whose 
action on monomials tp is defined by 

fxj) s * -^ s *(-i)Ei = i[('".+"'.)'+'».'»i]+9 (70) 
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As shown in appendix A, by using Eq. (66), it is not difficult to see that for 
Z we have also the representation 

Z = (dot A) J B cxp[-S*(^i)}d^ 1 ...d^ g , (71) 

where now the action S* is given by 

S*(ti) - l(lA~^) - V*(ti), (72) 

and the potential V* is related to V according to the following relationship 



f*(€,€) = iog 



where e = (1,1; 1,1). (73) 



Let 77,77 be vectors of Q g (Grassmann sources). If we define the generating 
Grassmann functional integral 

2(77,77) - J E 'exp[-S(€,|) + (7,,!) + {t,fjM 1 dS 1 ...dS B dS g , (74) 

we have also 5 

2(77, 77) = (dct A) J B eMS*(tl,ri,fiM 1 d^...d^ g , (75) 
where now S* is defined as before, but V* depends also on 77 and 77 and is given 

by 



V*(£,£;fj,ri) = log 



(76) 



Hence we now have, for example, that a fermionic system whose action 
is represented by a free part A plus a short range interaction, as regard the 
functional integral, can be viewed as a fermionic system with free part A^ 1 
plus a "long range" interaction. Substantially V* is the complement of V in 
Qig- However notice that the term "long range" here has not exactly the same 
meaning used in physics. The above long range character is related to the fact 
that if V connects few Grassmann variables, V* connects many. 

5 Let us observe that we adopted for Berezin integration the definition (3), which is different 
from J B ■ ■ ■ d£ g F(£) = fi k as ^> besides £1^2 • • -£g, depends on others Grassmann 

variables too. In this case it is important to consider the relations (2). 



24 



Acknowledgments 

I would like to thank G. Jona-Lasinio and C. Presilla for their encouragement 
and many useful discussions and K. Scharnhorst for critical comments on a 
previus version of this paper and for bringing Rcf. [11, 12, 10] to my attention. 
I would to thank also S. Tanase Nicola for communicating to me his derivation 
of eq. (12) and for interesting comments. 

APPENDIX A 

Let be 

a generic action. Since (£, A£) is an even element of Qi g , we have exp[— £)] = 
exp[— A£)] exp[V(£, £)] and we can expand the latter factor in Taylor series 
with respect to V^(£,£). Remembering that [V r (^,f)] Ar is finite sum of monomi- 
als, we see that Z can be written as a finite sum of Gaussian Berezin averages 
each one multiplied by a suitable coefficient given by the coefficients of the 
potential (68) and by detA (the Gaussian normalization). 
Explicitly Eq. (66) reads as 

(i\W(S) 

<4iti •••W .4 ' Si Si ■••? 9 9 ? 9 >A- X dctA ■ 

The r.h.s. of these relations tell us that we can calculate each Gaussian averages 
considering A -1 instead of A as matrix associated to the Gaussian action, pro- 
vide we use s* instead of s. ■0 s = £l~ ni £l~ ni ■ ■ ■ £g~" 9 £g~™ 9 can be obtained 
from ■0 s = i^" 1 ^™ 1 . . .£g 9 Cff 9 by multiple Grassmann derivation applied to the 
highest grade monomial ip e = £i£i-..<!; s ,£ 9 : 
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As regards (— l) w ( s \ it can be expressed in terms of s* = (mi, mi; . . . ; m g ,m g ) 
and one finds 



(— = (_i)W( s *( s )) = (_i)(^+™i)+2(m2+ro 2 )+...3(m 9 +m 9 ) (77) 
Hence, we have 

<^,0>A= E <*'>A% = E <^>^(-l)" (S * (S)) de7A 

SG(Z 2 *) 2 s SG(Z 2 *) 29 

and by using the definition of T it is easy to see that one has 

<na)> A =t<v { ^»'> A . 1 ^ J (78) 

and in general 

< V-(U) > A = f < > A -r (79) 

Finally by using Eq.s (78, 79) and the definition of V* (73) we get Eq. (71). 
By using this result the expression (75) follows as corollary. 

APPENDIX B 

Concerning (— l) Pw , since it accounts of the exchanges needed to bring in 
front of § 2 ^ 1+e *' ) and £j in front of £ ?( 1+ej ) and considering the structure of Cij, 
we see that the first permutation involves (g — 1) + (<? — i) exchanges and the 
second one g — j. We have 

(-l)fy = (-ljf-*-^- 1 . 

Thus for any given e e T 

S (e ij ) = A 2g (e ij )(-l)^- 1 A ij . 

Let us now calculate A 2g (e^ ) . Due to the arrangements of the generators of the 
algebra, we see from (16) that we can break up its expression in one part which 
contains only the barred components and another which contains both: 

A 2g(eij) = n(^r^ + ^T^ 1 ---^- 1 ) n(^T^ + i^ei . . . e fl ei . . . ej_i 
l=i ^ ' i=i^ 
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For the structure of e i3 we obtain 

A 2ff (6 ij ) = (-1)^ X (-1)3-1+9-; = { _ 1)3 -i-j-l 

Therefore we finally have 

s ( e ij) — Aij- 

APPENDIX C 

Among all the 2g factors which are contained in Ci g , only the i— th and the 
(.9 th can be different from 1, and one has 

ft,^,-(-i^)- 1 + ( -f"" + '-f"' 

x(-l)^ + --- + ^ 5+Wsl+ --- +w » _1+ ™ i+ --- + ™ 3+ ™ 1+ --- + " J - 1 (-l) £ ' +: ' +i+:L 

i.e. 

cwk, -c-d***) - 1+( -f + " 1+( -f' + "' (-ir^- 

x cxp{i7r[AT s l + ... + N% + Nl + ... + N 3 ^ 1 + n t + ... + n g + n x + ... + rij-i]}, 
where, due the structure of and by using Eqs. (31) and ( 32), one has 

1=1 
1=1 

APPENDIX D 

Let us suppose that = 1, let us consider in (49), the factor with labels (ij) 
and let us observe the structure of Hij. Let N[ lk = 1, we must consider four 
events: 

1. h<i + l, k>j — 1 from which none contribute to Hj from N hk 

2. h > i + 1, k > j - 1 from which in appears 9(s lj - s hk ) 

3. h < i + 1, k < j — 1 from which in Hj appears 0{s l i — s hk ) 
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4. h>i + l,k<j-l from which in appears 6(s lj - s hk ) + 8(s lj - s hk ). 

The case n. 1 gives (—1)° = 1. The case n. 2 involves (— iy( s ' J - shk )- : as one 
considers in (49) the factor with labels h, k one has i < h + 1 and j < k + 1 
from which j < k. But two processes cannot jump in the same row or the 
same column, it is effective only j < k — 1 and that means that in Hhk one 
has the analog of case n. 3 where i,j — > h,k and vicecersa. Then one has 
9{s hk — s*- 7 '), which together with #(s y — s hk ) coming from the former factor, 
give us (— l) 1 = —1. The case n. 3 is exactly analogous to the case n. 2. The 
case n. 4 gives (— l) 2 = 1. 

Therefore the integrals over the jump times are fictitious and it is then 
convenient to define the variables with four labels 9^ k satisfying the following 
properties: 

+ = 1 

Hence we will write: 

n* = [ E E w ; + E E 

ra—i+l l—l m— 1 l — l 

which does not depend on the jump times. 
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